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ABSTRACT 

In this paper we use the "Millennium Simulation" to re-examine the mass assem- 
bly history of dark matter halos and the age dependence of halo clustering. We use 
eight different definitions of halo formation times to characterize the different aspects 
of the assembly history of a dark matter halo. We find that these formation times have 
different dependence on halo mass. While some formation times characterize well the 
hierarchical nature of halo formation, in the sense that more massive halos have later 
formation, the trend is reversed for other definitions of the formation time. In partic- 
ular, the formation times that are likely to be related to the formation of galaxies in 
dark halos show strong trends of "down-sizing" , in that lower-mass halos form later. 
We also investigate how the correlation amplitude of dark matter halos depends on 
the different formation times. We find that this dependence is quite strong for some 
definitions of formation time but weak or absent for other definitions. In particular, 
the correlation amplitude of halos of a given mass is almost independent of their last 
major merger time. For the definitions that are expected to be more related to the 
formation of galaxies in dark halos, a significant assembly bias is found only for halos 
less massive than M*. We discuss our results in connection to the hierarchical as- 
sembly of dark matter halos, the "archaeological down-sizing" observed in the galaxy 
population, and the observed color-dependence of the clustering strength of galaxy 
groups and clusters. 
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1 INTRODUCTION 

The cold dark matter (CDM) scenario has been proven 
a very successful model for structure formation. In this 
framework, most mass in the cosmic density field ends 
up in virialized objects, called dark matter halos, and 
luminous objects, such as galaxies, are suppose d to form 
and evolve in such halos l|White fc Rees I Il978t) . Clearly, 
the first step in understanding the formation of galaxies 
is to understand the formation of the dark halo popula- 
tion. Using both TV-body simulations and semi-analytical 
methods, many important results have been obtained 
regarding the properties of dark matter halos in the 
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momentum p r opert y (e.g. [Barnes fc Efstathiou 
Cole fc Lacevl 1 19961 ; IBullock et al.l l2001al; IChen fc Jing I 

fe-g. iMo fc White I Il99' 



and clustering property 



Lemson fc K auffmann 1999; Shct h fc Torman 



1999tlsheth. Mo fc Torman 1120011 )! These results have been 



playing an important role in our understanding of galaxy 
formation. 

One important property of a dark halo is its forma- 
tion history. In the previous studies, this formation his- 
tory is usually characterized by a single parameter which 
is the time w hen a halo has acquired half of its final 
halo mass (e.g.lLacey fc Cole II 19931: Lemson & Kauffmann 



halo mass function (e.f 
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1993; Sheth & Torman 1 


19991; ISheth. Mo & Torman 112001; 


Warren et al.l |2006|), density and sub-halo profile (e.g. 


Navarro. Frenk & White 1 Il997l; IBullock et all 12001b; 



E-mail: liyun@astro.umass.edu 



1999 1; Ivan den Bosch 1120021 ; IGao et alj|2005l ; IWechsler et~ 
20061 ). This definition of halo formation time is useful be- 
cause it indicates when the main body of a halo is as- 
sembled. However, it is unclear if such definition is closely 
related to how gal a xies f orm i n a halo. F or ex ample, 
Ivan den Bosch et all (|2003l ) and lYang et all <|2003h both 
found that dark halos with masses around 10 11 ' 5 /i _1 M« 
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have the lowest mass-to-light ratio, which suggests that star 
formation is the most efficient in halos with a fixed mass 
around 10 n ~ 12 h' 1 M Q . Thus, for halos with masses much 
larger than this mass, the half-mass assembly time may have 
little to do with how galaxies may have formed in such halos. 
Based on the half mass formation time, more massive ha- 
los are expected to form later due to hierarchical clustering. 
This is in contrast with the recent observations that the stel- 
lar population in more massive system s are generally older 
|Thomas et al.ll2005l ; iNelan et al.[|2005h . This phenomenon, 
known as the "archaeological down-sizing", appears to be 
in contradiction with the "hierarchical" formation scenario, 
but may also indicate that the growth of galaxies in a halo 
does not follow the growth of the halo. 

Recently Gao, Springel fc White (2005, see also 
Wechsler et all 120061; lHarker et al.l [20061 : iJing, Suto fc Mo I 
20071 ; iGao fc Whiteil2007l ') found that, the half-mass assem- 
bly time of a halo is also correlated with halo clustering 
properties on large scales. Using iV-body simulations, these 
authors find that, for halos of a given present mass that is 
smaller than M*, the ones that assembled half of their fi- 
nal masses earlier are more strongly clustered in space. On 
the contrary, for halos more massive than M, , the ones that 
assembled half of their final masses later are more strongly 
clustered. If the star formation history is somehow correlated 
with dark halo formation history, as is expected from cur- 
rent theory of galaxy formation, these results would indicate 
that galaxy systems, such as clusters and groups, of the same 
mass but containing different galaxy populations should also 
show different clustering properties. Obs ervationally, there is 
evidence to support such connection (e.g. IWang et al . 2007b; 
I Yang et al.ll2006l ; lBerlind et al.ll200c3 ). Although there is still 
discrepancy among the different results, most studies show 
that redder groups are more stro ngly clustered than bluer 
groups. In particular, the results of lWang et al.l |2007bl ) sug- 
gest that such color dependence is stronger for groups with 
lower halo masses and becomes insignificant for groups with 
halo masses above ~ 10 13 /i _1 Mg. This dependence has the 
same trend as the assembly-time dependence of halo clus- 
tering, and it is tempting to link these two types of depen- 
dence. However, as mentioned above, the half-mass assembly 
time may not be a good indicator of the typical formation 
time of stars in a ha lo. Indeed, using a "shuffling" technique, 
ICroton et ail l|2007l ) found that, the dependence of halo bias 
on halo half-mass assembly time can only account for about 
half of the clustering bias seen in red halos in their semi- 
analytical catalogue. Clearly, in order to understand the ob- 
servational results in terms of halo assembly bias, one needs 
to define halo formation times that are more closely related 
to the formation of galaxies in dark matter halos. 

The main goal of this paper is to systematically study 
when various characteristic events take place in the halo as- 
sembly process and how they are correlated with halo mass 
and with the large-scale environments. To this end, we define 
a number of formation times to characterize each halo forma- 
tion history. We study in detail how each of these formation 
times is correlated with halo mass and how the halo correla- 
tion amplitude depends on these formation ti mes. Our anal- 
ysis is based on the "Millennium Simulation" (|Springel et al.l 
l2005al ). The paper is organized as follows. In Section 2, we 
briefly describe the simulation and the techniques to iden- 
tify halos and to construct merging trees. In Section 3, we 



describe our definitions of halo assembly times and how to 
estimate them from the simulation, and study how they are 
correlated with halo mass. In Section 4, we present the re- 
sults on the assembly-time (age) dependence of halo clus- 
tering. Finally, in Section 5, we summarize and discuss the 
implications of our results. 



2 THE SIMULATION 

In this paper we use the "Millennium Simulation" (MS) car - 
ried out by the Virgo Consortium (|Springel et alj l2005al ) . 
This simulation follows the evolution of 2160 3 dark mat- 
ter particles in a cubic box of 500 h~ 1 Mpc on a side. The 
particle mass is approximately 8.6 x 10 8 /["'M©, which en- 
ables us to study the assembly of halos more massive than 
around 10 11 /i _1 Mq with a reasonable mass resolution. The 
simulation adopts a flat ACDM with Qm = !ldm + = 
0.205+0.045 = 0.25, where fldm and Sib stand for the current 
densities of dark matter and baryons respectively; the linear 
r.m.s. density fluctuation in a sphere of an 8/i _1 Mpc radius, 
as, equals 0.9; and Hubble expansion parameter h = 0.73. 
In total there are 63 snapshot outputs between z = and 
z = 80, which are almost evenly placed in ln(l + z) space. In 
the MS simulation, the characteristic collapsing mass, M*, 
defined through cr(M*) = 1.69, is about 6 x 10 12 ft, _1 M Q . In 
order to identify dark halos, the Friends-Of-Friends (FOF) 
algorithm with a linking length b — 0.2 times the mean par- 
ticle separation is used, so that the structures identified (we 
will call them FOF groups hereafter) have a density approx- 
imately 200 times the mean cosmic density. In addition, by 
smoothing the FOF groups outside-in, each FOF group is 
also assigned a corresponding "virial halo" with a "virial 
mass" M v , so that the average density contrast between the 
"virial halo" and cosmic critical density p c , is 

A v (z) = 18tt 2 + 82[On(z) - 1] - 39[O m (z) - l] 2 (1) 

an fc Norman | [T998). The radius at which the density 
contrast first reaches A v (z) defines the virial radius, R v , of 
the halo. Despite the fact that M v is slightly (typically 5%) 
smaller than the corresponding FOF group mass, there are 
no other significant differences when one studies the accre- 
tion history of dark halos. In what follows, we always use 
"virial halos" in our analysis. Given a cosmological model, 
we define the virial velocity Kdr(z) of a growing halo at red- 
shift z as 



Kir(z) 



GM v (z) 
Rv(z) 



Av(z) 



1/6 



[M v (z)H(z) 



1/3 



(2) 



where R v (z) and H(z) are, respectively, the halo virial ra- 
dius and Hubble expansion parameter at redshift z. With 
this, we can follow the growth of the virial velocity (which is 
a measure of the halo gravitational depth) using the growth 
of the halo mass. 

The merging trees of dark halos in the MS are 
constructed on the basis of sub-halos. In each FOF 
group, self-bound sub-struc tures are identified using SUB- 
FIND (|Springel et al.ll200lf ). A sub-structure 1 at redshift 
21 is considered a progenitor of another sub-structure 2 at 
Z2 (zi > Z2) if a certain fraction of its most bound particles 
are in sub-structure 2. In addition to sub-halos, each FOF 
group contains one and only one "main virial halo", with 
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mass M v . In our analysis, we use M v to construct the mass 
growth history of the final halo. Since in general M v ac- 
counts for the central part of an FOF group, this treatment 
naturally avoids the ambiguous case where some accidently 
linked sub-halos that do not belong to the halo also con- 
tribute to the halo mass. 



3 HALO FORMATION TIMES 

As mentioned above, the formation of dark halos is a very 
complicated process. There are two ways to follow the mass 
growth of a halo with time. First, one can start with a halo 
at the present time, pick up the most massive progenitor 
in the adjacent snapshot at higher redshift, and repeat this 
procedure until the halo mass is so small that it cannot be 
resolved anymore. The mass accretion history of the halo is 
then represented by the growth of the progenitor mass along 
the "main branch". Alternatively, one can always look up 
the most massive progenitor at each redshift on the merging 
tree, and concatenate on these progenitors chronologically. 
The mass accretion history so obtained reflects the growth of 
the "maximum progenitor" . Note that the halo main branch 
does not always represent the maximum progenitor of a halo 
at each redshif t, especially at the early stage of a growing 
halo (e.g. iGao et all l2004bh . Most of the previous studies 



follows a power la w of M v (0), with power-index 7 = 0.32 (see 
lYang et al .11200 table 1). We therefore have, 



branch" when 


study- 


Lacev & Cole 


19931; 


Wechslcr et al. 


2006). 



ing the halo assembly histories (e.p 
Ivan den Bosch1l2002l ; lGao et alfeOOa 
In this paper, we use both definitions. We define the follow- 
ing set of parameters to characterize the assembly history of 
a halo: 

(i) £1/2, mb : This is the redshift at which the halo main 
branch has assembled half of its final mass, M v (0). This 
formation time has been widely used in the literature, as 
mentioned before. 

(ii) Zi/2 ti ; This is the highest redshift at which half of 
the final halo mass is contained in progenitors with masses 
(M p ) greater than 0.02 M V (0). The same kind of f ormat ion 
time has been used in iNavarro. Frenk fc White I (|l997h to 
characterize the formation time of a halo and to study how 
halo concentration is correlated with formation time. 

(iii) 2i/2,t 2 : This is the highest redshift at which half of 
the final halo mass has assembled into progenitors more 
m assive than a fixed m a ss, M c = 10 11 5 fe'^Q. A s show n 
by Ivan den Bosch et al.l (|2003l ) [see also lYang etafl i|2003h ]. 
halos with masses ~ M c have the minimum mass-to- light ra- 
tio, and thus are the most efficient in star formation. With 
M c = 10 11 ' 5 /i _1 M0, £i/2,t 2 therefore indicates when star 
formation starts to prevail the halo assembly history. By def- 
inition, only halos more massive than M c have a well-defined 
z l/2,t 2 - This formation ti me is analogou s to th e formation 
ti me, -ZN06, introduced by iNeistein et al.l l|2006h . According 
to lNeistein et all l|2006h . zno6 is the time when the sum of 
the progenitors above a given minimum mass, reaches half 
of the present day halo mass. 

(iv) Zm/l : This is the redshift when the progenitors more 
massive than M c have assembled a fraction / of M v (0) . Here 
the definition of / is based on the non-constan t mass-to- 
light ratio of dark matter halos (|Yang et al.l l2003). For halos 
more massive than M c , the mass-to-light ratio, M v (0)/L, 



f = a 



M v (0) 



Mv(0) 
Mc 



(3) 



where a is a constant, which is set so that / = § for M v (0) = 
M c - Thus, zm/l essentially reflects the time when a halo 
becomes capable of forming a fraction of its total stellar 
mass. Note again that zm/l can be defined only for M v (0) > 
Mc 

(v) z CO rc,mb: This is the highest redshift at which the 
halo's main branch assembles a mass of M c . This forma- 
tion time therefore indicates when a halo is able to host a 
relatively bright central galaxy. 

(vi) Zcorc.m P : This is the highest redshift at which the 
most massive progenitor has reached the mass M c . Note that 



for massive halos, z c 



may be different from 



(vii) z vmax : This is the redshift at which the halo's virial 
velocity V v i T (z) reaches its maximum value over the entire 
mass accretion history. According to equation ((2} , the value 
of Vvir(z) is expected to increase (decrease) with time, if the 
time scale for mass accretion is shorter (longer) than the 
time scale of the Hubble expansion. Therefore, z vmax indi- 
cates the time when the halo mass accretion tr ansits from 
a fast accretion pha se to a slow accretion phase l|Zhao et al.l 
l2003l ; lLi et"aill2007l ). 

(viii) zi mm : Last major merger time. Here we define a ma- 
jor merger as the event when the mass ratio between the 
smaller halo and the main halo is no less than 1/3. The last 
major merger time is defined to be the one when the last 
major merger occurred on the main branch of an assembling 
halo. 

Once the merging history of a halo is given, it is quite 
straightforward to determine the formation times defined 
above. The only exception is zi mm . Since the mass trans- 
fer from the merging halo to the main halo is a gradual 
process, a merger in general takes several snapshots to com- 
plete. Thus, if we used the halo mass increase in one time 
step, we would find only a small number of events in which 
the increase in the halo mass in a time step is large enough 
to be qualified as major mergers. In order to circumvent this 
problem, we start from one snapshot, and trace the progeni- 
tors (including those of sub-halos) back to all the snapshots 
within a 1-Gyr interval. As long as there is a progenitor 
with mass exceeding 1/3 of the main-branch halo mass at 
the same time, a major merger event is identified. The choice 
of 1 Gyr is not crucial; our tests using 0.5 Gyr or 1.5 Gyr 
give almost the same results. 

As illustration, we plot in Fig. [1] an actual merging his- 
tory of a typical halo selected from the MS simulation, with 
all formation redshifts defined above marked. As one can 
see, the different definitions give very different values of the 
formation redshift, and they capture quite different aspects 
of the assembly history of a dark matter halo. 

In Fig. [2] we show each of the formation redshifts as a 
function of halo mass. In each panel, the solid line represents 
the median in each mass bin, while the dashed lines repre- 
sent 20% and 80% percentiles, respectively. As one can see, 
less massive halos generally have higher values of z 1( / 2 , m b, 
z i/2,t 1 , z vmax and zi mm than massive ones, i.e. these for- 
mation redshifts have a negative correlation with halo mass. 
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Figure 1. Merging history of a typical MS halo, with all the defined formation times marked. Progenitors greater than 4 X 10 10 /i -1 M0 
are output at selected redshifts to avoid crowdness. The radius of each circle is roughly proportional to M 1 / 3 . Circles filled with hatch 
lines that are 45° clockwise to the vertical represent the main branch progenitors; while those filled with hatch lines that are 45° 
counter-clockwise to the vertical represent the maximum progenitors. 



Since these formation times are defined in a self-similar man- 
ner, i.e., do not involve any particular mass scale, it is not 
surprising that they show a similar "bottom-up" trend, a 
consequence of hierarchical clustering. Nevertheless, Zi/2,mbj 
z i/2,ti , Zvmax and zimm still represent quite different epoches 
of halo formation history, which can be seen from their dif- 
ferent values and scatter. For all halo masses, both z vmax and 
Z\mm have scatter that is much larger than Zi/2, m b, z i/2,t 1 - 
This indicates that both z, max and zi mm are more sensitive 
to the details of the halo assembly history. 

On the other hand, the other four formation redshifts, 
z i/2,t 2 , Zcorcmb and Zcorc.mp, all show positive correla- 
tion with the halo mass, in the sense that more massive halos 
experience these events earlier. For massive halos, z ore,mb is 
lower than z corc , mp and has larger scatter, which is due to the 
fact that for some massive halos, the most massive progen- 



itors are not in the main branch. The trend is particularly 
strong for z M/L , z cor e,mp and z cor c,mb- A halo of 10 12 /i _1 M Q 
assembles a progenitor of mass 10 n ' 5 /i _1 Mo typically at 
z ~ 1, while such a progenitor forms at z ~ 5 for halos with 
masses > 1O 14 /i _1 M0. Since z m /l, Zcorcmb and z cor c,mp are 
the redshifts that characterize when a halo was able to host 
a relatively bright galaxy, the results shown here suggest 
that massive galaxies can form much earlier in massive ha- 
los than in low-mass halos. If star formation in these mas- 
sive galaxies was eventually quenched as their stellar masses 
reach to some value, as is the case in the current AGN feed- 
back model, or as they merge into a mas sive halo where 
radiative cooling becom e s inefficient (e.g. IChurazov et al.l 
120051 ; ICroton et all l200fj| ; ICattaneo et alj |2008h . one would 
expect that the star formation activity shifts with the pas- 
sage of time from high-mass systems to the low-density field. 
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Figure 2. Formation time v.s. halo mass. Solid line represents the median of each mass bin, while dashed lines represent 20% and 80% 
quantiles, respectively. Note that in the panels for Zi/2,ti an d Zvmax, the slight drop at the low mass end (to the left of the vertical 
dashed line) is due to the finite mass resolution of the simulation, because small progenitors cannot all be resolved for halos that are too 
small. 



This may be related to the observed "down-sizing" effect 
that massive galaxies in present-day clusters in general have 
old stellar populations with little star formation activities, 
and most star formation activities at the present time have 
shifted to low-mass systems. This shift is perfectly consis- 
tent with the hierarchical formation of dark matter halos, 
provided that there are some mechanisms that can quench 
star formation in massive galaxies. As we have shown, more 
massive halos indeed assemble their masses later, but the 
formation of massive galaxies can actually start earlier in 
their progenitors. 

As mentioned before, the formation redshift zi/2,t 2 
defined h ere is similar to the formation time zno6 intro- 
duced by iNeistein et al.l (|2006l ) . However, the halo mass- 
dependence we obtain here is quantitatively different from 
theirs . At the massive end, the results of INeistein et al.l 
(2006) show continuous increase of the formation red- 
shift with halo mass, wh il e our s show a flattened relation. 
Note that INeistein et al.l l|2006l ) used the extended Press- 
Schechter formalism to generate halo merging trees, while 
we obtain halo merging trees directly from A r -body simula- 
tion. We suspect that the discrepancy may be due to the 
inaccuracy of the extended Press-Schechter formalism. 

In Fig. [3] we show z cor c,mb and z CO rc,m P versus zi/2, m b 
for halos of different masses. For low-mass halos, M v ~ 
lO 12 fr _1 M , z core ,mb and z corc , mp are very similar to Zi/ 2 ,mb- 
However, for halos more massive than 10 13 /i _1 Mq, z C ore,inb 
and Zcore.mp are both higher than 21/2, m b- in particular, for 
halos with M v > 1O 14 /i -1 M0, z C orc,m P ~ 7, without depend- 
ing strongly on the half-mass formation redshift, 21/2, m b- 
This shows again that, for massive halos, the progenitors 
that can host massive galaxies can form much earlier than 
when the halos assemble most of their masses. Thus, al- 



though dark halos form hierarchically, star formation may 
appear "anti-hierarchical" at late epochs when many halos 
in which star formation was efficient have merged into mas- 
sive systems. 



4 THE DEPENDENCE OF HALO 

CLUSTERING ON FORMATION TIMES 

Halos are biased tracer of the dark matter distribution. On 
large scales the auto-correlation function of dark halos, £hh, 
is expected to be parallel to that of the mass, £ mm , so that 
one can write 



£hh (M,r) = 6 2 £ mm (r), 



(4) 



where b is the so-called halo bias factor. Analytical mod- 
els and A r -body simulations have shown that the halo bias 
factor depends strongly on halo mass. Halos more massive 
than M, are more strongly clustered than the underlying 
mass (i.e. b > 1), while low-ma ss halos are less clustered 
than the mass (i.e. b < 1) (e.g. iMo fc White I [l996l: Ijine 
1 19981 ; ISheth fe Torman II 19991 ; ISheth. Mo fc Torman 1120011) 
More recently, using A r -body simulations. iGao et al.l 



found that for halos with a fixed mass, the halo bias fac- 
tor b also depends on the time when the halo first assem- 
bles half of its mass, i.e. on z 1 ^ 2 ,mb, in the sense that halos 
with higher Zi/2, m b are more strongly clustered (i.e. have 
a higher bias factor). This assembly bias is found to be 
stronger for halos of lower mass. Subsequent investigations 
using different simulations have confirmed this result (e.g . 



Wechsler et "all |2006|; lHarker et all 120061 ; IZhu et al.l 120061 ; 

Jing. Suto fc Moll2007i r and theoretical models have been 



proposed to understand the origin of such assembly bias (e.g. 
IWang. Mo fc Jing I l2007al ; ISandvik et aill2007l ; lHahn et all 
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Figure 3. Compare z COIeul ^ and z CO re,mp with Zi/2,mb- Each line represents the result for halos with the mass marked aside. For each 
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standard deviation. More massive halos generally have lower Zj /2,mb but higher z coro ,mb or z COIC>mp . Error bars in the right panel are 
generally smaller than in the left panel, indicating that the maximum progenitor could substantially deviate from the main branch, 
especially for massive halos at early time. 



l2007l : lDesiacaues 1120071 : iKeselman fc Nusser \\200l\ ). In most 
of these earlier investigations, the assembly bias is analyzed 
in terms of the half-mass assembly time, Zi/2,mb- However, 
as we discussed above, although zi/2, m b may be a good quan- 
tity to characterize the formation of the main body of a halo, 
it does not characterize other aspects of the halo formation 
histories that may be more closely related to the formation 
of galaxies in halos. With the various formation times we 
have obtained in this paper, it is interesting to investigate 
how the clustering of halos depends on these different for- 
mation times. 

In order to study the formation-time dependence of halo 
clustering, we use a 3-dimensional Fast- Fourier- Transform to 
derive the two-point correlation function of dark matter ha- 
los as well as dark matter particles on large scales. We then 
estimate the halo bias b for a given halo mass, using the 
square root of the ratio of the two-point correlation func- 
tion of halos and that of dark matter, averaged over data 
points in 8 h~ 1 Mpc < r < 20 fi _1 Mpc. This interval of r 
is chosen to ensure that the clustering is still in the quasi- 
linear regime where the linear bias relation Q is a good 
approximation. In Fig. [4] we show how halo bias depends on 



the formation times we have defined for halos of different 
masses. The dashed line in each panel shows the bias factor 
of the oldest 20% population among all halos as a function 
of halo mass, while the solid line shows the corresponding 
result of the youngest 20%. For comparison, we also show, 
as the grey line, the results for the total population with- 
out separation according to formation time. The errorbar on 
each data point is estimated through the error propagation 
function based on the Poisson noise of each data point of the 
halo two-point correlation in 8 /i -1 Mpc < r < 20 Mpc. 

As one can see, the bias factor of the total population 
increases with halo mass, and the increase is more rapid 
at the massive end. This result of mass-dependence of the 
halo bias factor is in good agreement with the results ob- 
tained in earlier investigations. In addition, for halos of a 
given mass, the bias factor also depends on the formation 
redshifts, although the strength of the dependence is not al- 
ways the same for different definitions. T he result based o n 
Zi/2, m b is very similar to that obtained bv lGao et al.l l|2005() . 
even though the result here is based on "virial" halos while 
Gao et al. used FOF groups. With the exception of the case 
with 2i mm , where no significant age dependence is found for 
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Figure 4. Age dependence of halo bias. Formation time used is indicated in each panel. Dashed lines are for oldest 20% halos while 
solid lines are for youngest 20% halos; the thick gray lines represent the bias of all the halos regardless of their ages. Error barshows the 
Poisson error. 



any halo masses, the strength of the age dependence in gen- 
eral decreases with increasing halo mass. For halos more 
massive than 10M, , we do not see any significant difference 
between halos of different ages. However, the noise here is 
too large due to the small number of systems available from 
the simulation, and so our data is not able to reveal the 
weak reversed trend, namely that the youngest halos are 
more strongly cluster ed, at the very massive end seen in 
some simulations (e.g. Ijing. Suto fc Mc~ll2007l ; IWetzel et alj 
120071 ). The strongest age dependence is seen in the cases 
°f £i/2,mb> z i/2, ti an d Zi/2,t 2 - Note that all these three for- 
mation redshifts are based on the properties of half of the 
final halo masses. This suggests that the assembly of the 
main parts of halos, especially the low-mass ones, may be 
affected significantly by large-scale environments. On the 
other hand, for the definitions that are based on the forma- 
tion of a progenitor of a fixed mass, such as 2 C oro,mb and 
Zcorcmp, the age dependence is weaker, particularly for ha- 
los with masses much higher than the progenitor mass, M c , 
used in the definition. As shown earlier, such progenitors 
in massive halos usually form at high redshifts where the 
large-scale environmental effects may not yet have time to 
develop. The age dependence based on z vmax is also weaker 
than that based on Z1/2, m b, presumably because the halo 
density involved in defining 2: vm ax is relatively high and so 
the the mass assembly before z vmal is less affected by the 
large-scale environment than that before zi/2, m b- Fig.|4]also 
shows that there is almost no dependence of the bias fac- 
tor on 2i mm . This suggests that major mergers may be con- 
trolled by the properties of the local density field, without 
being strongly modulated by lar ge-scale environments . This 
is consistent with the result of IPercival et al.l |2003), who 
found that, for halos with very recent major mergers (within 



the past 10 8 years), there is no difference in the halo bias 
compared with all the halos of similar mass. 



The formation-time dependence of the halo bias pre- 
sented above may have important implications. Previous 
studies suggest that the halo assembly bias may introduce 
obse rvable effects in the large-scale cluster ing of galaxies 
Ce.g. lNeistein et al.|[200rj| ; ICroton et al.ll2007t ). Using a large 
group catalogue constructed from the SDSS, IWang et al.l 
|2007br i found that groups with a redder central galaxy or 
a redder average c olor of member g alaxies show stronger 
clustering (see also lYang et a.1.1 |2006T > . Because of the com- 
plexity of halo assembly, it is unclear which aspects of the 
halo formation history are more closely related to the col- 
ors of the galaxies that form in halos. By "shuffling" galaxies 
contained in halos o f similar mass or formation time Zi/2, m b> 
ICroton et al. (2007) found, in their semi-analytical model, 
that the z^.mb-dependence of halo clustering can account 
at most half of the clustering bias of red galaxies. This im- 
plies that the difference in zi/2, m b alone may not be suf- 
ficient to account for the colors of galaxies. This result is 
not surprising, because the assembly history of a halo is 
quite complicated and it is not expected that 21/2, m b can 
provide a full characterization of such history. As demon- 
strated above, each of the eight formation times defined in 
this paper catches a different aspect of the halo formation 
history. It would be interesting to if see some combinations 
of these formation times are better correlated with the prop- 
erties of galaxies. As we have shown, the assembly bias be- 
comes insignificant for halos more massive than 10 13 /i _1 M o 
for most of the definitions of the assembly times. This may 
be the reason why the color-dependence of galaxy group 
clustering is si gnificant only for g roups less massive than 
~ 10 13 /i _1 M (|Wanget al.ll2007bl ). 
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Our results also show that there is virtually no de- 
pendence of halo bias on zi mm , the redshift of last major 
merger. In the literature, it has been suggested that ma- 
jor merger ma y effectively shut off the star formation in 
a galaxy (e.g., Hernquist II 19891: [Mihos fc Hernquistl 19961 ; 
ISpringel. Pi Matteo fc Hernquist Il2005bl ; iKang et al.ll2006h . 
and hence zi mm should be correlated with the current color 
of the central galaxy. However, if major mergers were the 
main reason to make a galaxy red, there would be no color- 
dependence of the clustering amplitu de of galaxy groups , 
contrary to the observational results of lWang et al.l (l2007bh . 
This suggests that major mergers alone cannot explain the 
red color of central galaxies. It is possible that the redden- 
ing of a central galaxy is the accumulative effect of a series 
of events triggered by, say, minor mergers, rather t han a 
dominant major merger fe.g. iGeorgakakis et al.ll2008l ). 



recent observations. However, since there is virtually no de- 
pendence of halo clustering on zi mm , the typical color of 
galaxies in a halo is not expected to be determined by the 
last major merger time of its host halo. 
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5 DISCUSSION AND CONCLUSIONS 

In this paper we examine the complexity of dark halo mass 
assembly history using the MS and using different formation 
times to characterize the various aspects of halo formation 
histories. We find that, formation times defined according to 
the assembly of a fixed fraction of the halo final mass char- 
acterize the hierarchical clustering, in the sense that halos 
of higher masses on average have later formation time. On 
the other hand, formation times defined by the formation of 
progenitors of a fixed mass where star formation is expected 
to be efficient, clearly show "anti-hierarchical" behavior, in 
the sense that halos of higher masses have earlier formation 
time. If some feedback processes can terminate the star for- 
mation in these progenitors, we would expect that galaxies 
in massive halos are redder, consistent with observation. We 
would also expect that the star-formation activities should 
shift with time from high-mass to low-mass halos, and so the 
observed "down-sizing" in star formation is not in conflict 
with hierarchical clustering. 

We also study how the clustering of dark matter ha- 
los depends on the various formation times defined in this 
paper. We find that halo bias shows a strong positive correla- 
tion with halo mass, in good agreement with earlier results. 
For fixed halo mass, our results confirm a positive correla- 
tion between halo f ormation time, Zi/2.mbi an d halo clus - 
tering strength (e.g. iGao et al.|[2005l ; IWechsler et~ai1l2006h . 
The strength of this dependence increases with decreasing 
halo mass. For halos more massive than 10 14 Ii~ 1 Mq, we do 
not find a clear reversal of the assembly bias. In general, 
for halos less massive than 10 14 /i _1 M0, there is a positive 
correlation between the various formation times defined and 
halo clustering strength, with the correlation being stronger 
for lower halo masses. However, the correlation amplitude is 
quite different when different formation time is considered. 
The strongest age dependence of halo clustering is seen on 
z i/2,t ± an d zi/2,t 2 - There is virtually no age dependence of 
halo clustering on halo last major merger time, zi mm , and 
the dependence on Zm/l, Zcorc.mb, z cor0imp is moderate. If 
the typical age of stars in a halo is correlated with halo as- 
sembly history in some way, then halos with fixed mass but 
containing redder member galaxies are expected to be more 
strongly clustered, and this color-dependence is expected to 
be weaker for more massive systems. This is consistent with 
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